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ABSTRACT 


The frequency response was determined for a second order servo- 
mechanism with backlash in the system. Frequency response for sinu- 
soidal inputs were obtained by digital computer solutions of the 
equations of motion of the system which had inertia and viscous 
friction on both sides of the gear train. System parameters included 
input amplitude, magnitude of backlash, damping ratio and resilience 
of the gear teeth. 

The Control Data Corporation 1604 Computer was utilized and 
Fortran programming was employed, 

Results show the effects of the amount of backlash and the 
gear material resilience on the output to input ratios, The exist- 
ence of sub-harmonic resonance was found with large values of back- 
lash. 

The authors express their gratitude to Dr. George J. Thaler 
of the Department of Electrical Engineering for his guidance and 
encouragement and to Professor Douglas G. Williams of the U. S. 
Naval Postgraduate School Computer Center for his generous assist- 


ance, 
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1. BACKGROUND AND INTRODUCTION 


Most servomechanisms have a mechanical connection between the 
driving components and the driven or positioning unit. Usually this 
connection is a gear train which, in most practical cases, has a 
finite amount of backlash in the gearing mechanism. This inherently 
renders the servo system non-linear and requires analysis techniques 
which are based on approximations o: which necessitate use of laborious 
graphical methods which may rely on simplifying assumptions. The 
backlash between the gear teeth cannot always be eliminated without 
possibly introducing other, perhaps less desirable effects such as 
coulomb friction and/or increased viscous friction. Manufacturing 
costs may also prohibit consideration of gears with mo backlasn. 

The relative distributions of viscous friction and inertia on 
each side of the pearing have direct effects on the behavior of the 
load, Graphical phase plane and analog simulation analysis techniques 
ао not reaqilv adapt to practical consideration and adjustment of 
these distributions. Also, they cannot handle the case of non- 
plastic impact of the gear teeth. This has led to the solution of 
this problem by the use of digital computers, 

Several investigations have been made of second order servo 
behavior using the digital computer at the U, S. Naval Postgraduate 
School. These include, limit cycle investigation by New for plastic 
impact (Ref. 1) and by Anderson and Luckett for plastic and elastic 
impact (Ref. 2), and a transient response investigation for plastic 


impact by Andrews and Kelley (Ref. 3). 


This thesis extends the coverage of this problem to a steady- 
state frequency response investigation, with both plastic and elastic 
impact, for a periodic (sinusoidal) input of varying frequency and 
amplitude. The general approach to the investigation was similar to 
that used by Anderson and Luckett. The equations of motion of the 
system were normalized so behavior would be independent of the un- 
damped natural frequency of the system, The parameters inertia, 
viscous friction, elastic "bounce" coefficient and damping ratio 
were in the same ranges of values used by Anderson and Luckett so 
that a correlation of results could be made if desired. 

Computations were made on a CDC 1604 digital computer using a 
Fortran System program which had definite advantages over the machine 
language programs used in the previous investigations mentioned above. 
Solution of the system equations was based on the phase plane approach 
as to whether the system was operating combined or separated. Even 
though the system was normalized, the computer program was made gen- 
eral to permit extension of the values of the parameters if desirable. 

This thesis includes investigation of the effects of the amount 
of backlash, the amplitude of the periodic input and the amount of 
system damping for both plastic and elastic impact of the gear teeth, 
The only non-linearity considered was backlash. 

The total system was piecewise linear in that one linear differen- 
tial equation described the motion of the system when the driving motor 
was combined with the load and two different linear differential eq- 
uations described the separated motions of the motor and load. 
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Boundary conditions for the differential equations were based on phase 
plane slopes and on equations based on the laws of conservation of 
energy and momentum which were related to whether gear contact was 
plastic or elastic. 

The original intention was to determine the effects of different 
distributions of inertia and viscous friction between the driving unit 
(motor) and the driven unit (load). However, because of the great 
amount of digital computer operating time required, it was necessary 
to confine the investigation to one ratio of inertias and one ratio of 
frictions. 

Results from the digital computer calculations provided data from 
which the magnitude of the output and the phase relation between input 
and output could both be determined, This thesis is concerned only 
with the magnitude of the output as compared with amplitude of the in- 
put sinusoid. | 

Special effects or phenomena such as jump resonance and sub- 
harmonic resonance, which are known to exist under certain conditions, 
were looked for in the course of the investigation. 

In the following sections the system equations are developed and 
the evolution of the computer program is described, A print of one 
version of the Fortran language computer program and the computer 


solution flow chart are included as appendices. 





2. SYSTEM DESCRIPTION AND DEVELOPMENT OF EQUATIONS 


2.1 Description of System 


For the purpose of this investigation, a second order system was 


chosem, The block diagram of the system is as follows: 


f ил n Ө 
SEK III E c 
1 2 1 


Motor f Load 
Gears 


Fig. 2-1 


This system is the same system that was reported by Anderson and 
Luckett in Ref. 2, Except for the non-linearity, the system is a 
simple second order servomechanism. There is only one non-linearity 
present, and that is backlash. The backlash is present in the gear 
train which separates the motor from the load. The effect of the 
backlash is to allow movement of the load and/or the motor inde- 
pendent of the movement of the other, so long as this movement does 
not exceed the amount of backlash. Also included as part of the gear 
train is the gear ratio; however, this has no special effect on the 
system other than is normally expected, 

Separated by the backlash is the motor inertia and friction and 
the load inertia and friction, The distribution of these quantities 


cannot be changed in the describing equations of the system due to 
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the backlash. In the linear case, the inertia and friction can be 
lumped on either side of the gear train without affecting the system 
operation, but this cannot be done when there is backlash present. 
The motion of the system can be described by a combination of 
three linear differential equations. When the gear teeth are touching, 
a single differential equation describes the motion. Once the gear 
teeth separate, one linear differential equation describes the motion 
of the motor, and another describes the motion of the load. Upon 
recombining of the gear teeth, the single differential equation again 
describes the motion. At the instant of separation or recombination, 
the appropriate differential equations have as initial conditions the 
last values of the equation which described the motion prior to the 
separation or recombination. The change over from one describing 
equation to the pair of equations and vice versa will be discussed 


further in the equation development section. 


2.2 Assumptions 


Prior to developing the equations which describe the motion of the 
system under consideration, the assumptions which have been made will 
be presented. These assumptions are essentially those made in Ref. 2 
and are repeated for convenience but with some modification. 


The following assumptions were made: 


1. The gear teeth were initially in contact, and the initial con- 
ditions of the system were all equal to zero. 


2. Plastic deformation of the gear teeth during steady state con- 
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tact or on impact and any torsional deformations of driving 
shafts are negligible. 

The inertias of the gears and drive shafts are considered as 
part of the load or motor inertia depending on their attach- 
ment in the system. 

The law of conservation of energy is completely satisfied in 
the perfect elastic case, and the law of conservation of 
momentum is satisfied in the imperfect elastic cases includ- 
ing the plastic case. When only the law of conservation of 
momentum is required to be satisfied, the energy lost by the 
system is dissipated in the form of heat by infinitesimal 
deformations of the gear teeth. 

The gear teeth are in contact only instantaneously during im- 
pact for all cases except the plastic case, and the impulse 
torques of drive, friction, bearing supports, etc., are zero 
during impact. 

The coefficient of restitution of the two opposing gear teeth 
is the same or can be described by an equivalent coefficient 
if the gear teeth are of unequal coefficients. 

Backlash is assumed to be equal at all points on the gear cir- 


cumference, and the backlash is measured at the output shaft. 


Equation Development 


From Fig. 2-1, the error signal, E, is the driving signal. This 


is equal to: 





Тһе torque developed by the motor, T, is equal to: 


“КЁ ze =$ Cu. (2-2) 


Equations (2-1) and (2-2) are true when the load and motor are combined 
and when they are separated. 


The load torque, Ti» 15: 


т. = 34,9 +f,0, (2-3) 


When the load and motor are combined, the load torque when referred to 


the motor side of the gear train is: 
Tim =P T, =P (8, +£,6,) (oe) 


The forward or positive direction of the system is defined to be 


when: 
O. =P Oy (2-5) 
and the backward or negative direction when: 
eG. = Powe (2-6) 


By differentiating either equation (2-5) or (2-6) with respect to 


time yields the same results. These are: 





SN s C (2-7) 


апа 
9. = の 6, (2-8) 


Substituting equations, (2-7) and (2-8), into equation (2-4) 


yields after rearranging: 


2 if = f e 
= 2-9) 
Im? Da tt д.) | 
The driving torque of the motor must equal the load torque referred to 


the motor side (combined system only) and the load torque of the motor 


itself. 


FE x ds Ir + е Shy + Tum (2-10) 


Substitution of equations (2-1), (2-7), (2-8), and (2-9) into (2-10) 


after simplifying gives: 


e o 
Sl Му а J КА | Bent Ge (2- №) 


(Un +2 J ) š T M 2 Ji " (Ju + の | 


Rewriting equation (2-11) 
= 2 2 
O. а түкте О + Wy O. = Why F (2-12) 


where 





EL 


Си = s 2-13 
n Jute J ( ) 


and 


fa tot, 
J Wh = — 
З Jv PR 


(2-14) 


Equation (2-12), therefore, defines the motion of the system when the 
gear teeth are combined, 

Whenever the gear teeth are separated, there is no driving force 
on the load, and its motion is defined by the following equation: 


"e f. 
oa J 


To be noted at this point is the change in notation for the load. 





O, = O (2-15) 


In order to distinguish between combined and separated operation, the 
subscript C or L is used respectively for the output. 
Combined operation: 


Output position - 


Sp 


Output velocity 


(|) 


(> 


Output acceleration - 
Separated operation: 
Output position - біз 


Qutput velocity 


CD: 
E 


Output acceleration 





5 
e 
j 
| 


Jn 





Although this convention in notation was adopted, it should be remembered 
that the same quantity is described by two different symbols, depending 
on whether the system is combined or not. 

By the gear teeth being separated, the load torque is not reflected 
back to the motor side of the gear train. Therefore, equation (2-10) 


reduces to: 


p š; i B, + ГА 8, (2-16) 


ог Бу substituting equation (2-2) into (2-16) and after rearranging 


reduces to: 


K 


м 


A 
Әм 


Since this investigation is directed towards the frequency response 


“ + . 7 
Um t 7 On + FO, = Ө, NEN 


of the system, the forcing function is: 


8. = À sin ot (2-18) 


If all angular frequencies are divided by the undamped natural angular 
frequency, W» this effectively is the same as setting w equal to 1, 
By this normalization, the angular frequency in equation (2-18) is 
really w/w . Therefore, the results are applicable to any simple second 
order system having any natural frequency Ws if the forcing function 
angular frequency is multiplied by the a of that system. 

Equation (2-12) describes the motion of the system when the gear 
teeth are combined. Equations (2-15) and (2-17) describe the motion of 


the load and motor respectively when the gear teeth are separated. 
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Left to be determined now is the instant when separation or combination 
occurs. 

Use of phase plane analysis is necessary to determine the instant 
of separation when the gear teeth are combined. By using velocity and 
position coordinates for the phase plane, the slope of a trajectory is: 

do 


lil = 18 (2-19) 





and specifically for the combined system: 


м Е. (2-20) 


_ а9 _ € 
№276 =, 


With proper manipulation of equation (2-12), the following obtains: 


2 2 ` 

М = on Qa - S G, 

= и 
Ө. 


(2-21) 


Using similar techniques to find the load alone slope on the phase 


trajectory from equation (2-15) yields: 





М = === (2-22) 
Ji 
Assuming that the system is combined, it will stay combined as 
long as the magnitude of the load alone trajectory slope is greater than 


the magnitude of the combined trajectory slope. Symbolically, the 


System will remain combined so long as the following is true: 
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> IN (2-23) 


С 


м 


Separation occurs when the combined trajectory slope equals the 





load alone trajectory. This, therefore, is the criterion to be satis- 
fied for transferring the motion description of the system from 
equation (2-12) to equations (2-15) and (2-17). 

For the system to recombine after separation, the gear teeth must 
come into contact with such velocities that they do not bounce apart in 
the elastic case. To recombine in the plastic case, the gear teeth 
merely must come into contact, 

In the general case, two criterion must be satisfied for the load 
to recombine with the motor, These are: 

l. The gear teeth must touch. 

2. The velocities of the load and motor must be such that they 

are equal after impact. (Elastic case only) 

Unless equation (2-5) or (2-6) is true, the gear teeth cannot be 
in contact. Therefore, when either equation (2-5) or (2-6) is satis- 
fied, the gear teeth have come into contact. In making this test, 
however, it must be remembered that O, is actually used in place of 

O. in equations (2-5) and (2-6). It is re-emphasized that these 
symbols actually represent the same quantity but are obtained by solving 
different equations. 

If the impact is plastic, the gear teeth will stay together. Ном- 


ever, in the elastic case, they do not necessarily have to remain in 


» 
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contact. During impact, the law of conservation of momentum must be 
satisfied; and after impact, the velocity ot the load and motor must 
be the same for the teeth to remain in contact. By using the lav of 
conservation of momentum, the respective velocities can be determined, 


In Ref. 2, the coefficient of restitution was defined. This is: 


_ (eó. - Š.) 
(278 
or rearranging: 


. / o Б 
= e O, = ۾‎ | 


ys 


where the primed terms represent the respective velocities after im- 


(2-24) 


pact. To satisfy the conservation of mementum law: 


EV 
ж 


must be true. То be noted, equation (2-25) is independent of the 


(9.-6,) = (6. - 6, (2-25) 


coefficient of restitution. 
・ / 
Solving equation (2-2h) for с” gives equation (2-2hA); and sub- 


stituting the result into equation (2-25), the following expression for 
* / 


= obtains: 


a = IR (PS. 


a a ко?) 28, (e+ | о. TN > (2-26) 
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If the gear teeth have come into contact, equations (2-26) and 


e / 
(2-24A) can be solved for c and Sr respectively. Unless 


Ak 


9, - 6 


М › the gear teeth will not remain together. 


In summary, assuming that the system starts in the positive direc- 


tion from zero with the gear teeth in contact, the operation is as 


follows: 


p 


Equation (2-12) defines the motion of the system until equation 
(2-23) is no longer true. 

Once equation (2-23) is no longer true, equations (2-15) and 
(2-17) describe the motion, taking as their initial values 

the terminal values of (2-12). In general, the direction of 
the system must be determined so that the proper value for 

©, can be used to initialize equation (2-17). ТЕ ). is 
positive, equation (2-5) applies. ТЕ O. is negative, 


equation (2-6) applies. In transferring from equation (2-12) 


e CASA 


During separated operation, recontact is determined from equa- 
tion (2-5) or (SO) 
In the event of recontact, equations (2-26) and (2-24A) are 
used to determine if bounce occurs. (Required only in elas- 
tic case) 
Finally, when there is no bounce, description of the system's 
motion reverts to equation (2-11). This time, the terminal 
value of Sm: used to initialize equation (2-12) along 

‚ / 
with Sg 
The steps of operation continue this pattern. 
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2.4 Parameters 

By referring to equations (2-12), (2-15), (2-17), and (2-24), the 
following parameters are available for varying: 

ee Т 

2. Wy 

5; er 

м 

5. fm 

EN 

Ье 

m e 

By setting GJ, equal to one, the system was made independent of 
its natural frequency. This. normalization was discussed previously. 
Not shown in the preceeding list as a parameter is the gear ratio, 
در‎ ‚ of the gear train. Since the gear ratio was considered to have 
no special effect on the system, its value was set equal to one. Also 
not shown is the magnitude of backlash. This was not used as a para- 
meter; however, it was assigned a value other than unity. 

The motor inertia and load inertia, Ju and Ji respectively, 
were combined such that with © equal to one their sum was one. 
(i.e., Ju + S Ji = 1) This means that any combination of J, and 


M 


л сап be used, зо long as their sum is one. 
By solving equation (2-13) for the system gain constant, K, the 


following is obtained. 


Wn LM +J.) 
yo 


15 


Kê 


(2-27) 


: : 2 7 
Since QO) is one, (Jy + p L9 is one and © is one; the gain con- 
stant, K, must also be one, This, therefore, was the value assigned to 
K in this investigation, 
Another variable was defined by combining the load and motor 


friction, f. and f 


L M? similar tọ the way the inertias were combined, 


This new variable called total friction is defined as follows: 


- 2 
= 2-28 
f, f, m ) 
The friction ratio of the system is defined to be EE 


; : | р В 
Substituting the value of one for w, / and (1, + の Jj) in 


equation (2-14), the following results 


ps = m A. (2-29) 


Or 


(2-294 ) 


tf 
5 


Assigning a value to the damping factor, J , fixes the sum of UE ES £, ). 
To determine the value of г, от Ж a choice also must be made for the 
friction ratio. Choice of Ü and the friction ratio, therefore, fixes 
fu and Е. 

For this investigation, the imput, ©. ; was a sinusoid defined by 
equation (2-18). The angular frequency of oscillation, w, was effec- 
tively normalized by making D equal to one. However, the values 
assigned to w were varied over the desired range of the frequency 


response. There is also no restriction on the choice of the forcing 


function amplitude A, 
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Equation (2-24) defines the coefficient of restitution of the s 
teeth. This definition imposes no restriction on its value; practi- 
cally however, it must be less than one and greater than zero. 

Having E many parameters available, it was necessary to restrict 
their variation. Listed below are the values assigned to the various 
parameters. It was felt that this choice would produce the most in- 


teresting results. The method used to vary parameters will be covered 


in Section 3. 


17 





PARAMETER 
Natural frequency 
Gear ratio 
Backlash 
Load inertia 
Motor inertia 
System gain constant 
Motor friction 
Load friction 
Friction ratio 
Damping ratio 
Forcing function frequency 
Forcing function amplitude 


Coefficient of restitution 


TABLE II-I 


Values Used 


for 


Parameters 


SYMBOL 
Wn 
P 
N 


L. 


Jm 


т 


oP ECL x 
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VALUES USED 


1.0 
1.0 
1,05 
.> 
.> 
1.0 
„12, „2, .36, 46, CO K. 
608, .16, .24, .32, 200 АС 
4 
105, 100.2, 9 о ое © 
ОТО Ó 
О 
QUE S 





3. DIGITAL COMPUTER PROGRAMMING AND METHOD OF SOLUTION 


ST Вир: Бат Computer Program 

А11 of the frequency response results were obtained by com- 
putations using the Control Data Corporation 1604 high-speed digital 
computer. The computer programming was done using the Fortran system 
with program input either by magnetic tape or by punched cards, Com- 
puter output was put onto magnetic tape and results were printed 
from the magnetic tape by use of an IBM 1401 high-speed line printer 
unit. 

The Fortran system was used because Fortran programming permitted 
inserting the equations of the servo system into the program in forms 
very close to their actual mathematical forms. The use of Fortran 
also allowed a great deal of flexibility in easily modifying the program, 
in changing system parameters, in inserting input data and in extract- 
ing the output results in a properly labeled decimal format. 

A print of a typical Fortran source program used in this inves- 
tigation is shown in Appendix B. This Fortran source program, includ- 
ing sub-routines, represents about 1600 machine language instructions. 

The system equations to be solved by the computer were developed 
in Section 2. The three differential equations describing the motion 
of the combined system and the motion of the separated motor and load 
were solved by subroutines which were parts of the overall Fortran 
program. These sub-routines were based on the Runge-Kutta numerical 
method for the solution of ordinary differential equations using a 


four stage process. Each ot the Runge-Kutta sub-routines utilized its 
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own integral sub-program. 

Even though the, system being investigated was normalized (a = 1) 
and a unity gear ratio was chosen いし の = 1), the program was generalized 
to provide for changing a and 2 to any desired magnitudes. Similarly, 
it was possible to select any value of system gain K, however a value of 
K = 1 was used for all calculation for reasons discussed in Section 2. 

As shown in Section 2, the forcing function or input to the system 
in this frequency response investigation was Ө, = A sin wt. The 
amplitude A of the sinusoidal input could readily be set prior to any 
computer run. The sinusoidal input was obtained from the Fortran 
Reference Library sine function which is an integral part of the Fortran 
system, The input frequencies w for each frequency response run were 
put into the program in an array of data. Each w was consecutively 
inserted into the calculations by the program until all desired ws 
were completed. A knowledge of Fortran programming is necessary to 
thoroughly understand how this is accomplished. 

It was necessary to choose integration intervals which would 
yield accurate results for the Runge-Kutta solutions of the differen- 
tial equations yet which would not require more calculations than 
necessary. It was found that for values of J 20.1 the following 


relationship was satisfactory for determining the intervals. 
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Integration interval At- . — (3-1) 
This resulted in 2000 x 2TY solutions per input cycle for all values 
of input frequency w. The values of At were also put into the pro- 
gram in an array of data and each value was inserted into the calcu- 
lation along with its corresponding w. 

With combination of low wand J <0.1, it was necessary to re- 
duce At to a value less than obtained from equation (3-1) in order to 
maintain sufficient accuracy, 

Since the damping ratio J for any run was dependent on the rela- 
tive magnitudes of > and f 


these parameters (`Y, Е, апа Е.) were 


L? M 


set prior to each separate run and were constant for that run. All 


runs were made with a constant ratio of 


pt, 


fu MES f. 


Even though the distribution of inertia between motor and load 
was even for all calculations, the program provided for changing this 
distribution to 80% of the total inertia to either the motor or load 
side. This was built into the program and arrays of values for Ju and 
Ji were put in as data similar to the arrays of w and At. 

The amount of backlash between the gear teeth was fixed before 
each computer run; in fact the value used for all but one run was 
A = 0.1 radian. A value of A= 0.05 was used in a special run to 
check the effect of changing A and А proportionally. 
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Because of the large number of solutions per input cycle (::00C TT ) 
and since it was not necessary to record all of these results, the 
program was designed to print the results for every one hundredth 
position computation. This provided sufficient results to obtain the 
frequency response of the system. 

Only the steady state response, after the system transients had 
effectively died out, was desired. Therefore a delay of five time con- 
stants (based on a similar linear system) was included in the program so 
outputs would not be ‘printed’ until this transient period had passed. 
The program mnenomic for this was DELAY. 

It was necessary to tell the program how long to compute before 
switching to the next value of w. This was done by causing the compu- 
tations to shift to the next w when a predetermined amount of problem 
time had elapsed after the transient period. This time was based on a 
desired number of input cycles. The program mnemonic for this was 


designated CYC. 


3.2 Sequence of Solution of System Equations 


The starting condition for each frequency response run was with 
the gear teeth in contact, i.e. the motor and load were combined, with 


this equation describing the system motion: 


O. + 27 o, O. + одг о. A sin wt (522!) 


This differential equation was solved as each increment of sinusoidal 


input was added. The corresponding value of Q. was calculated and the 
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combined phase plane slope was computed by 


ーー テー (3-3) 
ex 


This slope was compared with the computed value of the load-alone 


slope, which was constant for the run; 





N = (3-4) 


% 


When No = Ny the motor and the load were commencing to separate. If 


the system velocity ©, was positive at separation the gear teeth were 
separating from the original contact position, 


If O was negative, the teeth were separating from the opposite side of 


the backlash. 


After separation, the equatiou describing the motor motion was 


б, + 0, r e= A sin wt (3-5) 
M м M 


This differential equation was solved for each input increment. At 


the same problem time t, the following differential equation for the 


load alone was solved: 





9 + - 9 =O (3-6) 


From these solutions, the positions of the motor and load were 


compared in order to determine when the gear teeth had re-touched on 
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either side of the backlash, A condition to be satisfied for re- 
touching was 
OL = Om | (3-7) 
or 
ЕЕ (3-8) 


When the teeth did make contact, the load velocity after touching was 


% 


calculated from 


' | - fe, 
О, = ا‎ 9, (1+е) + 0, 176 (3-9) 


and the motor was 


ら = 9-e(66。-@) ғ” 


Ж 


The following equation was then solved to see if the teeth stayed 


combined 


Is J Op = O, LEY (3-19) 


Where € was an arbitrarily small number which was used to determine 
when the velocities were almost equal in the case of elastic impact. 
When the teeth stayed combined, the program returned to the 

combined differential equation (3-2) and the complete sequence was 


repeated. 


24 





3.3 Printing of Results . 

Solution of any of the three differential equations of the system 
produced the values of the respective position Ө апа velocity c 
Since the position and velocity of the combined system were defined as 
the load or output position and velocity, these values E oe. 
printed as the outputs when the gear teeth were ii. contact. When the 
motor and load were separated, the load alone was printed as the 
output, and the motor position and velocity were also printed for the 
same problem time t. 

For each print of system position and velocity, the problem time 
and the input position were also printed. Prints were made for each 
one hundredth solution of the differential equation or equations con- 
cerned. Special prints of problem time, position(s), velocity(s) and 
input were made just after the gear teeth separated or combined. 

A typical print of computer output is shown in Appendix C. This 
shows that output position and velocity (Ө апа Ө.) only were printed 
when the system was combined, while motor and load positions and 
velocities also were printed when the teeth were separated. 

Only five decimal places of results were printed because of 
limitations in the accuracy of plotting results. The computational 
accuracy was somewhat greater since the Runge-Rutta method was 
accurate to about (At), which was Io for the largest At used. 

A complete computer program flow chart and explanation is in- 


cluded in Appendix A, 
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3.4 Computation Time | 

Computer operating time to obtain a steady state output varied 
considerably for different values of J and for the different input 
frequencies w. For example, to compute through a five time constant 
transient period, with a low damping ratio of Ў = 0.1 (long time con- 
stant) and with an w of 2.0 (small integration interval), required 
about 28 minutes of computation. On top of this, it took about one 
and three-quarter minutes to compute and print the outputs for each 
input cycle. From one to eight cycles of input (after the estimated 
transient period) were required to obtain a steady state peak output. 
For the other extreme of a high X of 0.8 and a low a of 0.1, the 
total computation and print time required to obtain an output peak 
took about two minutes. 

It is apparent that the investigation required a great amount 
of computer time. All computer runs, including test runs, took a 
total of about 83 hours of computer operation. 

An attempt was made to develop a computer program which would 
approximate the energy stored in the system for steady state at the 
first input frequency and would retain this "stored energy when の 
was changed. It was hoped this would materially reduce the transient 
period in reaching steady state at the new у. Test runs were promising 
but time did not permit refinement of this program to where it could 


be used in the investigation. 
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4. RESULTS AND DISCUSSION 


4.1 Determining M from the Computer Output 
From standard frequency response nomenclature, the symbol M is 
defined to be the magnitude of the ratio of system output to system 


input. Symbolically this is: 


Ө. 


М = 6. 





(el) 


where 
Ө =А 


Appendix C shows a sample print out of data from the computer. 
O: 


it was necessary to pick off the value of output position at its maximum 








One of the quantities printed is output position. To determine Р 
or minimum. In general, there were two peaks. One was a positive peak 
and the other a negative peak. These two peaks came about because the 
input was permitted to go through only one cycle once printing began. 
The magnitude of these peaks were not always identical. This indicated 
that the system had not reached steady state; however, in most cases, 
the difference between the peaks was less than plotting accuracy. Тһе 
Ө. 
With values of the damping ratio greater than about .4, the method 


Ө, 


about O.3 and lower, there began to appear an oscillation at the natural 








average of the magnitude of these peaks was taken to be 


of obtaining as described above was used. At damping ratios of 








frequency of the system super-imposed on the output wave form. Its 


e 


ҮШ 





presence could be expected when a limit cycle would be predicted by 
Ref. (2). 
Fig. 4-la shows this phenomenon. The damping ratio, E: » was 

0.1 іп this case. Shown in Fig. 4-lb is the output waveform of the 
same system with a higher damping ratio. ( f 20.5) It can be 
seen here that there is no super-imposed "limit сус1е". To obtain 

le. for the system whose output. is shown in Fig. 4-la, it was nec- 
essary to take the mean of the super-imposed oscillation. In the case 
shown, this was not particularly difficult. However, as w was in- 
creased towards the undamped natural frequency, this procedure became 
more difficult. Fig. 4-2a shows the output waveform for the same 
system at = .6. From this figure, it can be seen that it is vir- 
tually impossible to establish the mean of the "limit cycle" oscilla- 
tions. Even if the output wave shape were available for a longer time 
period, it would still be virtually impossible to obtain an accurate 
mean value. То compare the output at the same forcing frequency, but 
with higher damping ratio, Fig. 4-2b was drawn. Here Р =0.5, and it 
can be seen that there is no difficulty picking the peak output value, 

Since the method of picking Ө varied, it was not possible to 

program the selection of | ©. | . In some cases, a programmed method 
would have been satisfactory, but unless all the data were examined, 
e 


the value of thus obtained would be questionable. 
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4.2 Frequency Response Variation with Damping Ratio 
Frequency response curves, for plastic impact of the gear teeth, 


А 


were obtained for values of Kof 10, 5 and 0.95, with damping ratio 

jf varying from 0.1 to 0.8. These curves are shown in Figures 4-3, 
4-4 and 4-5 respectively. For comparison purposes, a similar family 
of curves for a comparable linear system is shown in Fig. 4-6. 

It was determined by test runs, where both A and Á were varied 
proportionally, that the frequency response was identical as long as 
the ratio А маз constant, For example, the same frequency response 
resulted with A = 1.0 and A = 0.1 as with A = 2.0 and A = 0.2. This 
makes the frequency response results generally applicable to all 
simple second order servo systems for the same relative amounts of 
backlash since the system was normalized with œ = l. For this reason 
the abscissae of the curves are labeled w/w , and the curves can be 
applied to any simple second order system by multiplying by w of the 
system. 

The resonant frequencies, for values of damping ratio that resulted 
in resonance Е 0.6), меге found to increase as the amount of back- 
lash was increased. This increase in resonant frequency was greatest 
for the higher values of Y while it was almost insignificant at a 
ША ОС 0.1. 

The peak values of the ratio of output position to input ampli- 
tude at resonance мз increased above the linear peaks as the amount 
of backlash was increased. The increase was greatest at the lowest 


value of J$. A plot showing this variation from the linear for dif- 


ferent amounts of backlash is shown in Fig. 4-7. 
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As might be expected, the effect of backlash is greatest at the 
resonance frequencies. The frequency response curves, for all amounts 
of backlash, are above the linear curve at values of w below the re- 
sonant frequencies. At frequencies greater than the undamped natural 
frequency and for values of Y > 0.4, the response with backlash drops 
off more rapidly than the linear and the degree of drop-off is greater 
as the backlash is increased. In the same frequency range, the response 
curve at low values of Y is slightly above the linear curve (ignoring 
sub-harmonic resonance which will be discussed later). This indicates 
that with heavy damping, backlash results in a more rapid attenuation 
of output at frequencies above the undamped natural frequency of the 
system, while with light damping, attenuation is slightly less than the 


linear up to 20. 
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1.3 Effect of Magnitude of Backlash on Frequency Response for Constant E と 


A more detailed examination of the effect of the amount of backlash 
in the system compared with the input amplitude was made for a constant 
damping ratio. The value of J = 0.5 was chosen since it represents a 
practical damping ratio which is less than optimum and results in a 
moderately oscillatory system. Also, no limit cycle exists at this 
value of Y according to Anderson and Luckett, Ref. 2. It was pointed 
out earlier that as long as the ratio of A/A is constant the frequency 
response is invariant. Therefore, A/A is a logical system parameter. 

Figure 4-8 shows the frequency responses at 了 = 0,5 for A/A of 
10, 5 апа 0.95 with plastic impact of the gear teeth. The linear re- 
sponse is also shown on the plot. It can be seen that the response 
deviated from the linear even when the magnitude of backlash was only 
one-tenth of the input amplitude, The deviation increased as the rela- 
tive magnitude of backlash was increased. The effect of backlash was 
to change the response curve to resemble the response for a lower damp- 
ing ratio. 

At frequencies above the undamped natural frequency, with 
A/A = 0.95, the response curve dropped below the linear curve and 
reached a minimum near w= 2. Up to this input frequency the ratio 
of output frequency to input frequency was one to one. With a slightly 
higher w, the ratio of output to input frequencies f rapidly trans- 
sitioned to 1/2 and the magnitude of the output abruptly increased. 
This was the first of a series of sub-harmonic resonance peaks which 


were encountered as J) was increased further. The actual values of M 
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between sub-harmonic resonance peaks were not all definitive as will Бе 
shown in the discussion of the output wave shapes for this region of w. 
АЕ some frequencies where input frequency was not an integer multiple 
of the output frequency, the output magnitude could not be defined 
because of the random, unsteady nature of the load position. This is 


shown by a dashed line in Fig. h-8, 


4.4 Sub-harmonic Resonance with Plastic Impact 


The first sub-harmonic resonance peak occurred at w = 2.2, (Fig. 4-8). 
At input frequencies of 2.0 and below, the outputs were steady with ampli- 
tudes decreasing with increasing w, and the ratios of output to input 
frequencies were jue The wave shapes of the output for (9 - 1,8, 1.9 
and 2.0 are shown in Figs. 4-9 and 4-10. At w = 2.1 the output was 
very unsteady (Fig. 4-10b) and a value for M could not be determined. 

At this input frequency, the ratio of output frequency to input frequency 
a, was not constant but fluctuated with an average of between 1/1 and 
2. 

When w was increased to 2.2, the frequency ratio an became a 
steady 1/2. The output wave shape became steady but unsymmetrical, as 
shown in Fig. 4-lla. The peaks of this output wave occurred at the times 
of contact of the gear teeth. The reason for the unsymmetrical wave 
shape was due to the difference in relative velocities of the motor and 
load at the time of contact. Referring to Fig. 4-lla, the higher 
(positive) peaks occurred as a result of the lower relative velocity 
at the instant prior to the preceding contact (smaller peaks). This 


smaller relative velocity imparted a higher velocity to the load which 
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caused a greater excursion, The opposite was true in going from contact 
at the higher peaks (high la velocity) to the smaller peaks, This 
phenomenon is believed to be related to the fact that there was plastic 
impact of the gear teeth. The value of M at the sub-harmonic resonant 
peak at w= 2.2 in Fig. 4-8 came from the average of the high and low 
peaks. 

When W was increased to 2.3 the output wave became erratic, the 
ratio mt. was transitioning from 1/2 to 1/3 and the exact value of M 
was indeterminant but was obviously less then at t - 2.2 (Fig. 4-11b). 

At w= 2.4 the output wave (Fig. 4-12a) was similar to the wave 
shape at the previous sub-harmonic resonance peak of w= 2.2. The wave 
was steady but unsymmetrical and the ratio of an was a fixed 1/34 
The value of M was further reduced (averaged from peaks of output wave). 

Going to w= 2,5 (Fig. 4-12b) the wave shape was again unsteady 
but A was a consistent 1/3 and M began to increase. This increase 
in M continued for the next four ws. 

Wave shapes for ays 2.6, 2.7, 2.8 and 2.9 shown in Figs. 4-13 and 
4-14, show very convincingly a sub-harmonic resonance. At w= 2.6, the 
general wave shape of the output during this resonance is first observed. 
For all four of these frequencies, the system output had the same wave 
Shape. As the frequency was increased, the M value also increased. 

All four of these frequencies gave a n of 1/3. Just as the output 
wave shape was very steady, the p was equally as steady at 43. 
The sub-harmonic rise associated with these frequencies is shown on 
Fig. 4-8. Since the value of M for these frequencies was so well 


defined, this portion of the curve was not dashed as was required for 
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frequencies between w = 2.0 and w = 2.4. 

A further increase in frequency beyond w = 2.9 resulted in a drop- 
off in M, Fig. 4-15 shows the output wave shape for w= 3.0. This 
wave shape wastaken for a period of over 50 time constants. As can be 
seen, there are portions of this wave which are quite steady and 
portions which are very unsteady. The steady portion had a ove of 
Wie The M value taken from this steady portion was used to plot on 
Fig. 4-8. No M value could be obtained from the unsteady portion due 
to its inconsistent movement. The ЕЕЕ. during the unsteady portions 
was not defined either. It averaged between 1/3 to 1/h. 

At frequencies greater than 3.0, the output motion was erratic. 
This can be seen in Fig. 4-16 for w = 3.1 and 4.0 as well as in 
Fig. 4-17 for w= 5.0. Since M could not be deduced from these wave 
shapes, no points are plotted on Fig. 4-8 beyond w = 3.0. 

The БЕН for w = 3.1, Fig. 4-16a, was 1/3. This may be mislead- 
ing because only a small amount of the output wave was obtained, At 
ш = 4.0, Fig. h-16b, the ЕЕ was between 1/h and 1/5. Fig. 4-17, 
for = 5.0, shows the output wave shape for a considerable time. It 
can be seen here that the wave was not consistent. The РЕ: for this 
w varied widely. Depending on the portion of the curve examined, the 
ЕТЕ. was as small as 1/6 and as large as 1/4. 

Although no sub-harmonic resonance was found beyond w = 2.9, it 
would appear that it does exist. In transitioning from one sub- 
harmonic frequency to another, the output wave was quite unsteady. 
When a definite output frequency was established, the M was also 


defined. Since the AR ratio changed from 1/3 at w= 2.9 to 1/6 at 
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w= 5.0, a sub-harmonic resonance at 1/4 and 1/5 must have occurred between 
these two frequencies. 

From the limited data available, it appears that the system tends 
to resonate at odd sub-harmonics. Frequencies between W = 2.0 апа 
w= 2.2 did not have a fixed БЕСІ Бы Then at Q) = 2.2, a steady 1/2 ratio 
was observed. 

During the build up to the third sub-harmonic resonance peak, a 
steady ratio of 1/3 was present. This steady EER over several values 
of frequency, leads one to believe the system would tend to seek an odd 
sub-harmonic frequency. Immediately past the sub-harmonic peaks, the 
system began seeking the next sub-harmonic frequency. The even sub- 
harmonics did not appear steady except at w = 2.2. There was a con- 
siderable range, w= 2.4 to w= 2.9, where Ai was a steady 1/3 add- 


ing support to the belief that the system would seek odd sub-harmonics. 
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1.5 Effect of Elastic Impacu 

Frequency response data were obtained, at a constant Y = 0.5, for 
values of A/A equal to 10, 5 and 0,95 with e = 0.6 and e = 0.8. These 
parameters were the same as for the results discussed in Section u 
except for the values of the coefficient of restitution e, The response 
curves for e = 0.6 and e = 0.8 are shown in Figs. 4-18 and 4-19. The 
values of M for all three values of A/A were found to be slightly 
less for both cases of elastic impact than for e = O (Fig. 4-8) at 
frequencies near resonance. This indicates that elastic impact tends 
to reduce M as compared with plastic impact. There was practically no 
difference between the response curves for e = 0.6 and e = 0.8 up to 
input frequencies of W= 2. 

With А/А = 0.95 sub-harmonic resonance began to occur at about 
ш = 2.8 for e = 0.6. Up to w= 2,7 (Fig. 4-20a) the ratio EN was 
a steady 1/1 and the output wave shapes were steady and symmetrical, 
At w = 2.0 the output waveshape became erratic (Fig. 4-2la) and De nn 
was transitioning from 1/1 со o The response curve remained above 
the linear at all frequencies (Fig. 4-18), unlike the equivalent curve 
for e = O (Fig. 4-8) which dropped below the linear between 1. and the 
occurrence of sub-harmonic resonance. 

The response curve reached a minimum at w= 2.8. At the next 
test frequency W = 2.9, the output increased and the ratio К was 
between 1/2 and 1/3. No even harmonic was found as there was for e = 0. 
The wave shape at w = 2.9 was also uneven with very inconsistent peaks 


(Fig. 4-20b). At w= 3.0 and w = 3.2, ER was still between 1/2 and 
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1/3 the value of M was increasing with unsteady waveshapes (Figs. 4-2lb 
and 11-220 AE Ш = 3.8 (Fig. h-22b) ЕТІН became a consistent 1/3 and 
M remained about the same, but the wave shape was still not steady. As 
w was increased to 3,7 and to 4.0 (Fig. 4-23) the waves became symmetrical, 
BE. remained 1/3 and the value of M was decreasing. At the highest 
test frequency w = 4.3, M was still decreasing, the wave shape was be- 
coming unsymmetrical (Fig. 4-24) and “уг, was decreasing toward 1/4. 

In comparing the response curves of Fig. 4-8 (e = O and Fig. 4-18 
(e = 0.6) for A/A = 0,95, it can be seen that the addition of bounce 
changed the curve and the sub-harmonic resonance pattern at frequencies 
above ш. Bounce eliminated the drop ot the response curve below the 
linear response and apparently al o eliminated the even sub-harmonic 
resonance peak. The odd (1/3) sub-harmonic resonance peak was lower 
with bounce but was wider, i.e. the value of ics = 1/3 persisted over 
a wider range of frequencies. Sub-harmonic resonance started at a 
higher frequency when bounce was introduced. 

Еот е - 0.8 апа A/A = 0.95, sub-harmonic resonance started to 
occur at about the same value of w as for e = 0.6. The investigation 
did not go beyond this frequency for e = 0.8, 

In general, elastic impact of the gear teeth tended to reduce the 
effects of the nonlinearity for both the magnitude of the output and the 


amount of sub-harmonic resonance, 


4.6 Response Characteristics with Very Low Damping Ratio 


In an effort to find jump resonance and the factors influencing 


its occurrence, some runs were made with a damping ratio, Y , equal to 
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0.05. Considerable computer time was required for calculating the out- 
put at such a low damping ratio. Therefore, only a limited number of 
frequencies were investigated. 

The basic parameters for this investigation were: damping ratio 
equal to 0.05, A/A equal to 1.2, coefficient of restitution equal to 
zero, The frequencies used were between 0.5 and 1.8 with emphasis on 
the frequencies near the undamped natural frequency. Results of these 
runs did not show any jump resonance; however, some interesting things 
were shown. 

For a linear system, the resonant frequency is always less than the 
undamped natural frequency. (Resonant frequency is the frequency at 
which the resonant peak occurs on a frequency response.) Regardless 
of how small the damping ratio gets, the resonant frequency is always 
less than the undamped natural frequency; however, as Í approaches 
zero the resonant frequency approaches but never exceeds the undamped 
natural frequency. 

With backlash present, this is not always so. Fig, 4-25 shows a 
resonant peak at a frequency between 1.01 and 1.02 which was greater 
than (Do It has been shown previously that the effect of increasing 
backlash was similar to reducing the damping ratio. In addition to 
this effect, it appears that with a sufficiently large backlash and 
small enough J, the resonant frequency is increased more than the in- 
crease associated with a reduction in (ў alone, This may not be 
limited to low damping ratios. It may have occurred at higher damping 
ratios but was masked due to insufficient data being taken in the 
vicinity of the resonant peak. 
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The height of the resonant peak for this low Y is considerably 
higher than for the linear case. However, extrapolation of Fig. 4-7 
gives approximately the same value for this peak as shown on Fig. 4-25. 
At frequencies near the resonant frequency, between O,95 and 1.05, the 
output wave shapes were steady. The frequency ratio was also steady at 
one to one, Fig. 4+26 shows the output wave shape at w = 0.95 and 
Fig. 4-27 shows it at w - 1.0, Аз «ап Бе seen, these are steady. 

The output below œ = 0.9 was rather difficult to interpret. Ву 
the previously discussed methods for obtaining M, it was not possible 
to determine this quantity accurately. Therefore, values of M for these 
frequencies were not plotted on Fig. 4-25, 

Fig. 4-28 shows the wave shape for w= 1.2. Even though the output 
wave shape is steady at this frequency, the input frequency is not an 
even multiple of the output frequency, This uneven relationship has 
been mentioned previously (Section 4.4) as an indication of the onset 
of sub-«harmonic resonance. Therefore, the low damping ratio decreases 
the frequency at which sub-harmonic resonance appears, Insufficient data 
are available to determine the exact location of the first sub-harmonic 
resonant peak, but it appears to occur at a frequency near w = 1.6. The 
M at this peak is approximately 9, and the frequency ratio is even at 


one to two. 
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4.7 Test of Program and Validity of Results 


For there to be complete confidence in the results of this investi- 
gation, some mention should be made of the steps taken to prove the com- 
puter program from which the results were obtained. Also, some mention 
of the validity of the results is in order. It is the purpose of this 
section to discuss these subjects. 

Since there was no easy way that the system investigated could be 
simulated on an analog computer; and since little information is avail- 
able on such non-linear systems, it was necessary to test the program 
by approximation. First of all, the program was divided into two basic 
parts. These were, the combined system portion and the separated system 
portion. As integral parts of these two sections were the basic decisions 
which had to be made to transfer from one to the other. 

It was rather simple to test the combined section. This was done 
by setting the backlash, А » equal to zero, Doing this effectively 
eliminated the separated system portion and the decisions required to 
transfer into that portion. With the backlash equal to zero, the output 
should be the same as the output from a linear system. Comparison of 
the M obtained from the print out data, as described in section 4-1, 
with linear results obtained by conventional means were very favorable. 
The slight differences could be attributed to the computer solution 
still being in a transient state. The data from the computer solution 
was taken after five time constants had elapsed. The following table 


shows this comparison. 
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TABLE 1У-1 


Me 8 Ажо 


の M, computer | TT 
„l .9971 ル 4 .99716 
ES „98810 „98821 
= 97172 Оо 
n , ui 693 . 9 69% 
> .91180 _ .91192 


It was not possible to test the separated portion alone.  There- 
fore, it was necessary to test this portion as part of the overall test 
of the program. For a linear system, the ratio of forcing function am- 
plitude to the amount of backlash, A/A › 15 infinite since the backlash 
is zero, Therefore, if this infinite ratio is approximated by a non- 
linear system, the results should be very close to linear results. To 
approximate the linear system, an A/A of 100 was chosen. The results 
of this test run were within 2.5 percent of the linear response at all 
frequencies less than the resonance frequency. At frequencies above 
resonance, the non-linear response was less than five percent different 
than the linear response out to the highest frequency used. The maximum 
percent difference occurred at the highest frequencies where small dif- 
ferences account for higher percentages. 

Since the program was capable of approximating a linear system 
while actually going through all of the calculations necessary for a non- 


linear system; it was reasonable to assume that with a greater amount 
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of backlash, smaller the program would give valid results. 


=, 

Even though tte time solution of the differential equations was 
considered accurate, the method of determining M may have introduced 
inaccuracies, For damping ratios greater than 0.3, the method of 
obtaining M was quite accurate, The only improvement would have been 
to calculate for a longer time, but this would have required more time 
than could be justified by the increased accuracy. 

At damping ratios less than 0.3, the interpretation of the output 
waveform, る was not so clear cut. АЕ these damping ratios, the 
"limit cycle" was super-imposed on the output wave. At some frequencies, 
w less than about 0.4, the method used to obtain M was quite good. 

(Mean of "limit cycle" oscillations ) However, for frequencies between 
about 0.4 and 0.7, the M was not so nicely defined. It was in this 
frequency range for low 76 that the results were not as accurate ав 
they were otherwise. At all frequencies above an (p of about 0.7, the 
results wer« quite good except as dicussed in the sub-harmonic section. 

Mention has been made that data was taken after five time constants 
based on a linear system, The wave shapes at earlier times were looked 
at to determine if this delay were required. From these, it was evident 
that the system was in a transient state, Even after five time constants, 
there was evidence of transients. However, the transient affect was 
quite small. Depending upon the parameters, the transient was virtually 
gone shortly after five time constants had elapsed. The only exception 
to this was discussed as part of sub-harmonic resonance, In the cases 
where steady state was reached, it was reached in less than six time 


constauts. 
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4,8 Application of Results and Extension of Method 


while there is limited direct application of the results of this 
investigation to servo Jesign, there are some points that could prove 
helpful in the design of certain types of systems. If a system with 
backlash were to be operated with a periodic input, Fig. 4-7 would be 
useful in determining how much backlash could be tolerated for a given 
input amplitude to meet specified peak outputs, Also, the results 
showing the effect of gear resilience can be of benefit in choosing 
a gear material, since higher resilience tends to reduce the effect of 
the backlash. 

The general form of the computer program can be adapted to calcu- 
late any particular frequency response of similar second order systems 
J 


with known system constants (K, J A, etc.) for specified para- 


М? `1? 


meters. The program can be very easily modified to obtain transient 
response of second order systems with any desired sets of constants and 
parameters. Similarly, it can be readily applíed to steady state (limit 
cycle) ínvestigations. 

As has been pointed out, the method of obtaining frequency response 
results in this investigation requires a great deal of digital computer 
time. If further investigations were to be based on the same general 
approach, especially at very low damping ratios, a technique of going 
from one frequency to another without having to go through a new tran- 
Sient period each time, should be devised. Such a program was briefly 
attempted, as mentioned in Section 3, which transferred boundary con- 


ditions of the differential equations at the first frequency to the 
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appropriate differential equations as initial conditions at the next 
frequency. This was an attempt to utilize the "energy stored in the 
system to more rapidly set up а steady state condition for the next 
frequency. It appears this concept may be the only way that would make 
feasible this digital method to more extensive investigations. It could 
result in showing jump resonance when it exists, where the program used 


herein apparently did not. 
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5, CONCLUSIONS 


Based on the findings of this investigetion, the following conclu- 
sions are made: 

1. The system frequency response remains unchanged if A and А 
are varied proportionally, with all other parameters constant. 

2, The presence of backlash in the system has an effect similar 
to lowering the damping ratio of the system, 

3. The effect of backlash is more pronounced at lower damping 
ratios. 

4. Increasing the backlash causes a more oscillatory response. 

5. With a very low damping ratio and a small value of A/A » the 
first resonant peak can occur at a frequency greater than the undamped 
natural frequency. 

6. At sub-harmonic resonance the ratio of input frequency to output 
frequency is an integer. 

7. The output magnitude is erratic when the input frequency is 
not an integer multiple of the output frequency. 

8. The odd sub-harmonics are more pronounced; and with elastic 
impact, the even sub-harmonic does not appear. 

9. The lower the damping ratio, the lower the frequency at which 
sub-harmonic resonance will occur. 

10. With all other parameters constant, sub-harmonic resonance 
will occur at a lower frequency with plastic impact than with elastic 


impact. 


3 





11, The effect of elastic impact is to reduce the effect of back- 
lash. 

12. With a coefficient of restitution of either 0.6 or 0,8, the 
frequency response of the system is essentially the same, 

13, Combinations of low damping ratio and high frequency require 
excessive computer time to obtain the frequency response. 

14. When a limit cycle exists (Ref. 2), an oscillation at the 
undamped natural frequency of the system is superimposed on the output 


motion, 
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APPENDIX A - FLOW CHART OF COMPUTER OPERATION 








yes 


No 


yes 4 7 


48 





e No 


19 





TO. 


jll. 


l, 


14, 


19% 


16. 


EXFLANA¡ ION OF FLOW CHART NUMBERS 


Set initial conditions and assign values to coustants and parameters 
of system equations: еа f, s f , A, A e K, JO , WwW» 
Problem time tor delay of five time constants, number of input 
cycles desired after five time constants, 

XPrint heading for results. 

*Print run number and values of e, J and P 

21 column headings for Time, Ө,, С” Ө,» Ө, and Ss 


Set first initial conditions for combined differential equation: 


Time = О, E = O, Ө = 0 
C с 
Select values of «p and At. 


*Print values of Wy, At, and J,. 
L. 


Solve differential equation for combined operation 
um . 2 2 : 
0.--2У/шө,0, “ол O, = wn A sin wt 


Increment time by integration imntcrval ЕЕ + At. 
Calculate value of input 9 。 

Does time equal five time coustants? 

Is backlash zero? 


Solve for ,, 


e А 2 . 
6, =-2Yw,0, - 0,0, + wm Asınwt 
Calculate ^ 
N= ê, 


Calculate { 
L 


Does 
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кү. 
18. 


19. 
20, 


21, 


eT. 


28, 


29. 


50, 


EL. 


Count a time increment At. 

Gone through 100 at's’ 

Set At count to zero, 

*Print values of Time, Oc Ө. апа p 
Does time equal five time constants plus time for desired number 

of input cycles after transient period? 

Set initial conditions for second and third differentíal equations. 
Does time equal five time constants? 

*Print values of Time, Ө, 9 and Өр. 


Are gear teeth separating from original contact position plus A? 


Calculate Ө, = a 
9. 


Calculate oy = 





Wa 


Solve differential equation for motor alone 


Om K К | 
Br Ры 

J J 

M M M 
Solve differential equation for load alone 

О +50 =0 

ыы. 

Increment time t— + Ac. 
Calculate value oi input Ө. 
Does time equai five time constants? 
Count a time increment At, 
Gone through 100 at's? 
Set At count to zero. 
*Pri f Time , O0, O, and 
Print values of Time, Op 9, N? M an т 
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39. 


HO. 


hl, 


ho, 
ТЕЙ 
hh, 
h5. 
hé. 
87. 
48. 


* 


Does time equal five time constants plus time for desired number 
of input cycles after transient period? 

Are gear teeth touching in initial contact position? 

Are gear teeth touching in initial position plus A? 


ы 
Calculate load bounce velocity or 


5? “м 


S J. 


рг 901+ Ө те 


‚ / 
Calculate motor bounce velocity Ө 


e 


Ó = Eu CP Ou e on 
M fP 

Do gear teeth stay together? 
Set initial conditions for combined differential equation, 
Does time equal five time constants? 
*Print values of Time, 9,» 9 = Ө, апа z 
Go to combined differential equation. 
Completed calculation for all values of œw? 


Stop. 


| fi s と Е と 
Print refers to putting the indicated items or their values onto 
magnetic tape for later printing, 
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APPENDIX D 
COMPUTER FORTRAN SOURCE PROGRAM 
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